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Abstract. This paper reports the detection of a heavy element enriched RV Tauri variable with an abundance pattern that differs 
significantly from a standard s-process enriched object: V453 Oph. Based on optical high-resolution spectra, we determined 
that this object of low intrinsic metallicity ([Fe/H]= -2.2) has a mild, but significant, enrichment ([s/Fe] — 1-0.5) of heavy 
elements for which the distribution points to slow neutron capture nucleosynthesis. This result is strengthened by a comparative 
analysis to the non-enriched RV Tauri star DS Aqr ([s/Fe]= 0.0). Although V453 Oph is the first RV Tauri star showing a strong 
s-process signature, it is not accompanied by C enhancement, challenging our current nucleosynthetic models of post-AGB 
stars that predict a simultaneous enrichment in C and s-process elements. The low N abundance excludes CN cycling as being 
responsible for the low C abundance. We explore three different scenarios to explain the heavy element distribution in this 
evolved object: an enrichment of the parental cloud, an accretion scenario in which the chemical patterns were acquired by 
mass transfer in a binary system and an intrinsic enrichment by dredge-up. 
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1. Introduction 

RV Tauri stars are rare supergiants that occupy the high- 
luminosity end of the pop. II Cepheid instability strip. The 
identification of RV Tauri stars is based on their characteris- 
tic light curves, which show alternating deep and shallow min- 
ima with periods ranging from 30 to 150 days. This variability 
is often accompanied with lower amplitude irregularities and 
cycle-to-cycle va riability (e.g. review on pop. II Cepheids by 
IWallersteirJ2002l 

The far-IR excess of RV Tauri st ars due to circumstellar 
dust was detected by IRAS and led ljural dl986l) to conclude 
that they are young post-AGB stars, in which the circumstellar 
dust is a relic of the strong dusty mass loss on the AGB. For 
many RV Tauri var iables, a near-IR excess is detected as well 
(Gehrz 1972; Evans 1985). The genuine high-luminosity nature 
of the RV Tauri class was confirmed by the detection of some 
members in Globular Clusters and in the LMC. The LMC vari- 
ables obey a period-luminosity (P-L) relation, which follows 
the high luminosity extension of the P-L relation of W Vir stars 
jAlcocket al.ll998l) . 
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Euler telescope also on La Silla 



Chemically, however, RV Tauri stars do not show evidence 
for post-third dredge-up enrichment: no high C-abundances ei- 
ther s-process overabundances due to a third dredge-up have 
been found. A summary of the available quantitat ive abun- 
dance data on RV Tauri stars can be found in iGiridhar et alJ 
(2000J and lMaasetaD J2005I) . Many RV Tauri stars show ev- 
idence of selective depletion in which the refractory elements 
are depleted relative to the non-refractory ones. These photo- 
spheric chemical patterns are the result of a poorly understood 
process in which the circumstellar dust is separated from the 
circumstellar gas , followed by a selective re-accretion of th e 
non-refractories (Mat his&Lamerslll992t IWaters et alj fl992). 
Important diagnostics for this process are high [Zn/Fe], [Zn/Ti] 
and [S/a] ratios, because Zn and S have a very low condensa- 
tion temperature with respect to Fe, Ti and the c-elements . The 
process will also result in a C/O ratio larger than solar due to the 
stability of the CO molecule, which is thought to be abundant 
in the gas accreted onto the stellar surface. Since s-process ele- 
ments are refractory, a subsequent depletion process can mask 
an intrinsic AGB enrichment. Nevertheless, overabundances of 
s-process elements relative to other species with similar con- 
densation temperature were undetected. 

The depletion process can be very efficient in post-AGB bi- 
nary stars for which there is evidence that the dust is napped 

1 II "3 

in a circumstellar disc (Van Winckel et al. 1995). This geom- 
etry may well be a necessary condition for the depletion pro- 
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Table 1. Some basic parameters of V453 Oph and DS Aqr. 
The pulsation periods , spectr al types and magnitude ranges are 
from Kholop ov et alJ l[l998). The other parameters were ob- 
tained from Simbad. 



Data 


V453 Oph 


DS Aqr 


Other name: 


BD-02 4345 


HD 216457 


Coordinates: 


a-2000 +17 26 49.13 


+22 53 17.04 




£2000 -02 23 36.44 


-18 35 30.99 


Min. mag.: 


10.4 mag 


11.08 mag 


Max. mag.: 


11.53 mag 


11.90 mag 


Spectral type: 


Fp 


F2II 


Period: 


81.3d 


78.213 d 


Galactic 


1 20.72 


44.10 


coordinates: 


b 17.57 


-61.56 



cess to be efficient JWaters et al.lll992l) . Evidence is growing 
that also in dusty RVTauri stars the dust is trap ped in a sta- 
ble disc lEvanslll985t I Van Winckel et alJll999ah and the bi- 
nary nature of several dusty RV Tauri stars has been established 
(Van Winckel 2003, and references therein). It is, however, not 
generally accepted that all dusty RV Tauri stars reside in bina- 
ries. There is no clear relation between the current IR-exces s 
and the strength of the depletion process (Giridhar et al. 2000). 

The tradition al spectros copic classification of RVTauri 
stars bv lPreston et alJ fT*963') is strongly determined by the ef- 
ficiency of the depletion process. Spectroscopic group A con- 
sists of stars of spectral type G or K, showing strong absorption 
lines but normal CN or CH bands, while some also show TiO 
bands at minimum light. Group B stars, being generally some- 
what hotter, show weak metal lines combined with enhanced 
CN and CH bands. Members of group C resemble those of 
group B, except that enhanced bands of CN are not seen. It 
is now clear that the low metallicity of the weak-lined group B 
stars is determined by the efficiency of the depletion process. 
Depletion patterns are less strong, if present at all, in group A 
stars. Group C stars are not affected by depletion, probably due 
to their low intrinsic metallicity (Giridh ar et all20 00). 

In this paper we focus on the chemical analysis of two 
RV Tauri stars, V453 Oph and DS Aqr of spectroscopic class 
C to study potential AGB enrichment. The luminosity of these 
stars is about 2400 ± 1600 L , as determin ed using the P-L re - 
lation of RVTauri variables in the LMC jAlcock et alj fl998). 
This high luminosity identifies them as likely post-AGB stars 
(e.g. lBloeckerll995l) . Since class C objects are not affected by 
the depletion process, these intrinsically metal-poor objects are 
ideal to study AGB enrichment and possibly the metallicity de- 
pendence of the s-process. Some basic data of both stars are 
shown in Tabled 

After an overview of the observations, we discuss our spec- 
troscopic parameter determination in Sect. [3] The abundance 
pattern for V453 Oph and DS Aqr is derived in Sect.|4]and dis- 
cussed in Sect. [5] In Sects. 6 and 7 we explore different sce- 
narios to explain the observed abundances. We end with the 
conclusions in Sect.[H] 



2. Observations 

2.1. Spectroscopic 

High-resolution spectra of both RV Tauri variables were ac- 
quired with the fibre fed spectrograph FEROS, mounted on the 
1 .5 m ESO telescope at La Silla. The resolving power of the 
spectra is R ~ 48000 at 5000 A and a wavelength coverage 
ranging from 3600 A to 9200 A is obtained in a single expo- 
sure. Two exposures were obtained for V453 Oph and one for 
DS Aqr. The log of the observations is reported in Table|2] 

Table 2. The log of the obtained FEROS spectra together with 
the signal-to-noise (S/N) ratio in different spectral regions. The 
S/N ratio listed for V453 0ph corresponds to the final spec- 
trum. 

Name Date and UT S/N S/N S/N 

400 nm 585 nm 770 nm 



V453 Oph 28/06/01 05:30 
28/06/01 06:25 
DS Aqr 27/06/01 09:05 



50 
50 



180 
110 



180 
80 



The reduction of the spectra was performed in the standard 
fashion, using the FEROS DRS package in the MIDAS envi- 
ronment. Our reduction includes bias correction, background 
subtraction, cosmic hit cleaning, flatfielding, order extraction 
and wavelength calibration. Ultimately, the spectra of the dif- 
ferent orders were merged and normalised. 

For V453 Oph, the two normalised spectra were averaged, 
weighted by their signal-to-noise (S/N) ratio. In Table[5]the S/N 
ratio is given for different spectral regions. Sample spectra are 
shown in Figs.[5][6]and[7] 

Additionally, we monitored DSAqr with the CORALIE 
spectrograph mou nted on the Swiss telescope at La Silla, 
Chile. CORALIE JOueloz etalJll999h is a high resolution 
(R ~ 50000) fibre-fed echelle spectrograph specifically con- 
structed to obtain high precision radial velocities through 
a cross-correlation technique. We started the monitoring in 
August 2000 and obtained 40 data points. 

2.2. Photometric 



Using the light curve ephemeredes of iKholopov et alJ £998), 
we estimate that the spectra of V453 Oph and DS Aqr were 
taken at pulsation phase <p — 0.4 and <p = 0.0 respectively. 
Unfortunately, no simultaneous photometric data over a large 
wavelength range could be obtained for these phases. Instead 
we used photometric data obtained around the phases <p = 0.55 
and <p = 0.3 respectively. These data are reported in Table|3 

The total reddening was determined minimising the differ- 
ence between model atmospheres and dereddened fluxes in the 
optical wavelength range. For the model atmospheres, Kurucz 
models with parameters determined in Sect. |3]were used. This 
is possible because small photometric phase discrepancies are 
only expected to result in small shifts in atmospheric parame- 
ters. 
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For V453 Oph, we obtain a total reddening of E(B-V)= 
0.50 and for DSAqr, E(B-V)= 0.05. The resulting SED for 
both stars is shown in Fig.^ No near-infrared excess was de- 
tected for either variable. The high reddening of V453 Oph is 
consistent with the E(B-V) obtained from the diffuse interstel- 
lar bands at /l=5780 A and 5797 A, of w hich the equiva lent 
width correlates well with this parameter (Josafatsson & Snowl 
Il987h . Using the data reported in iKrelowski et alJ dl999t) 
these bands give an E(B-V)=0.48 and 0.54 respectively for 
V453 0ph. 



Table 3. The photometric measurements of V453 0ph and 
DS Aqr. The data are taken from various a uthors and plotted 
in Fig. G] a Data from lGoldsmith et alJ Il987h : b Data from the 
2MASS catalogue; c Data from lGehrzTdl972l) . 



band V453 Oph (mag) DS Aqr (mag) 



u 


11.76 a 




11.66 a 




B 


11.33 a 




11.18 a 




V 


10.54 a 




10.51 a 




R 


9.99 a 




10.12 a 




I 


9.39 a 




9.73 a 




J 


8.58 a 


8.806 b 


9.16 a 


9.488 b 


2.2 fim 




7.6 C 






H 


8.09 a 


8.442 b 


8.79 a 


9.197" 


K 


7.95 a 


8.283 b 


8.72 a 


9.089 b 


L 


> 7.0 a 




> 7.0 a 




M 


> 6.3" 




> 6.3 a 




N 


>4.5 a 




>4.5 a 




1 1.3 (im 




> 4.6 C 







3. Parameter determination 

3.1. Atomic data 

For the abundance analysis, all lines showing a clear sym- 
metric profile were measured. These lines were identified, 
based on wavel ength, using mainly the solar line list of 
iThevenir] dl989h . This list was also employed to identify 
blended lines using the calculated equivalent width (EW) of 
all lines in the spectral region. To calculate the abundances 
of the different elements, atomic data were preferably taken 
from the list with accurate oscillator strengths for A to F 
type stars, w hich has been collected at the Instituut voor 
Sterrenkunde (Ivan Winckel & Revniersl l200ol iRevniers et alJ 
2004). Alternatively, to obtain accurate oscillator strengths, 
the line lists of the Vienna Atomic Line Database (VALD2) 
were u sed (http://www.astro.univie.ac.at/~vald/ iKupka et"al] 
Jl999l) ). Additionally, the line list of the D.R.E.A.M. 
project (http://www.umh.ac.be/~astro/dream.htmll was used 
for Ce, while for La, Eu and Dy, data rep orted i n re- 
spectively lLawler et all lT2001al) . lLawler et aiT j^QOlb) and 
iBiemont & LoweHl993l) were used. 




V453 Oph 



10 2 10 3 10* 

wavelength A(nm) 
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10-9 
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Fig. 1. The SED for V453 Oph (top) and DS Aqr (bottom), con- 
structed using a total reddening of E(B-V)= 0.50, 0.05 respec- 
tively. The full line is the Kurucz model, while the other sym- 
bols are photometric measurements from various authors, listed 
in Table Upper limits are marked with arrows down. The 
plus signs are taken fro m the IUE archive, the diamonds from 
iGoldsmith etail dl987h . th e circle s from the 2MASS archive 
and the squares are from Gehr3 dl972l) . The figure shows 
clearly that no near-IR excess is observed. 

3.2. Radial velocities 

The heliocentric radial velocities for V453 Oph and DS Aqr 
were determined using Fe lines as -123.8 + 0.6kms _1 and 
-23.9 + 0.8kms _I respectively. These values differ respec- 
tively only by 5 km s~ l and 3 km s _1 from published values 
dGiridharetalJl998l). 

For DS Aqr, we have also 40 high resolution, low signal- 
to-noise CORALIE spectra (see Sect. 12. 1> . In the online reduc- 
tion of these spectra, radial velocities are determined by cross- 
correlation with a standard F0 mask. Since many of the corre- 
lation profiles were quite noisy, we decided to construct a mask 
optimized for DS Aqr. This technique is proven to be very ef- 
ficient for other RV Tauri stars dMaasl2003l) . The special mask 
contains only the lines that were used in the abundance study 
(Sect.0}. In Fig.|2]we compare the correlation profiles for two 
spectra obtained by using the standard mask and by using our 
constructed mask. Their smoother profile and their deeper cen- 
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i standard F0 mask 
■ DS Aqr mask 

^ u i , i i , I i 

-60 -50 -40 -30 
v rai (kms- 1 ) 



0.7 



-60 -50 -40 -30 -20 

V rad ( klrl S ~') 



Fig. 2. Cross-correlation profiles for DSAqr. Comparison of 
the cross-correlation profiles for two CORALIE spectra, ob- 
tained with a standard mask (open squares) and a specific tem- 
plate for DS Aqr (full squares). The profiles on the left panel 
are for the spectrum taken on 21/07/2001, the ones on the right 
panel for the spectrum of 25/07/2001. 

tral depth clearly indicate the higher quality of the latter. Six 
out of the 40 spectra were of a too low signal-to-noise to obtain 
a reliable velocity. 

The remaining 34 velocities for DS Aqr display a large 
peak-to-peak amplitude of 26 km s (see Fig. [3} Such a high 
value is not uncommon even for single RV Tau stars, and does 
not necessarily suggest a binary nature of the object. The photo- 
metric phase distribution of our 34 velocities is poor and since 
cycle-to-cycle variability is often observed in RV Tauri stars, 
we have as yet no significant test to probe the eventual binary 
nature of DS Aqr. 

Concerning the radial velocities of the two programme 
stars, we can conclude that the variability of the radial veloci- 
ties can be accounted for by the known photometric variability 
of the stars, with peak-to-peak amplitudes in V of 1.1 mag for 
V453 Oph and 0.8 mag for DS Aqr. 

3.3. Determination of atmospheric parameters 

The abundance pattern was com puted using the late st ATLAS9 
LTE model atmospheres (Castel li & Kuruca 120041) in combi- 
nation with the LTE line analysis program MOOG of Sneden 
(version 2002). The model atmospheres are uniquely deter- 
mined by the effective temperature (Trff), gravity (log(g)) and 
metallicity ([Fe/H]). They are computed with a constant micro- 
turbulent velocity (£, = 2 km s~') for the opacity distribution 
functions. 

We derived the photospheric model parameters solely using 
spectroscopic criteria. A temperature is determined by forcing 
the abundance derived from Fe I lines to be independent of their 
excitation potential (EP). The large range in EP and the very 
large number of lines (V453 0ph: EP = 0.0 - 4.4 eV, N FeI = 
104, N Fe n = 32; DSAqr: EP = 0.0 - 5.0eV, N FeI = 267, 
NFeii = 49) ensures an accurate identification of this tempera- 



DS Aqr 

CORALIE velocities 
folded on P = 78.2 days 



Fig. 3. Phase diagram of the 34 CORALIE velocities (filled 
circles) and the velocity obtained from the FEROS spec- 
trum (open circle) of DS Aqr. The velocities are folded on 
the photometric period of 78.213 days that was published by 
Kho lopov et al.l Jl998). Since cycle-to-cycle intrinsic variabil- 
ity is often observed in RV Tauri stars, a much better sampling 
is needed to be able to test for binary motion. 

ture. This is shown in Fig.|4]and confirmed by other elements 
with a large number of detected lines (e.g. Ti and Cr). 



0.2 



Cl. 

o 

» 0.1 



V453 Oph 
DS Aqr 




5000 



5500 



6000 
T ff (K) 



6500 



7000 



Fig. 4. The slope of the trend between the abundance derived 
from Fel lines and the EP versus temperature. The values for 
V453 Oph are depicted using squares and those for DS Aqr us- 
ing diamonds. These values were determined using log(g) = 
1.5, and log(g) = 0.5 respectively. An estimate on the errors 
was made by varying the micro turbulent velocity by ± 1 km s _1 . 

The surface gravity was determined demanding ionisation 
balance between Fe I and Fe II lines. The microturbulent veloc- 
ity was estimated requiring the independence of the abundance 
derived from Fe I lines of the reduced equivalent width. The re- 
sulting stellar parameters are shown in Table [4] where the Fe 
abundance is determined using an input model with a metal- 
licity of [M/H] = -2.0 and [M/H] = -1.5 for V453 0ph and 
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DS Aqr respectively. We used the conservative errors on the 
model atmosphere parameters of AT e ff = 250 K and A log(g) = 
0.5 (see Fig. @}. The error on the microturbulent velocity is 
adopted to be Atf, = 0.5 km s _1 . 

Table 4. The stellar parameters of the stars, as determined us- 
ing only spectroscopic criteria. The reported error on the Fe 
abundance is only due to line-to-line scatter. See text for more 
details. 

Name T eft (K) log(g) ^(kms' 1 ) [Fe/H] 

V453 Oph 6250 O lO -2.23 ±0.12 

DSAqr 5750 0.5 3.5 -1.62 ±0.12 



4. Abundance analysis 

4.1. Homogeneous analysis 

Severe selection criteria were used to incorporate or discard an 
absorption line in the analysis. The selected absorption lines 
have a small equivalent width (EW< 150 mA) and a clear sym- 
metric profile. Blended lines were only taken into account if 
clearly present and unblended in the other star. The abundances 
derived from these lines are shown in Table [5] 

The abundance of the CNO elements is always difficult to 
determine, particularly for such low metallicity stars. For C, the 
well-known multiplet at 71 15 A is undetectable for both stars. 
Therefore, a spectrum synthesis was performed using four C 
lines around 9070 A of which the result is shown in Fig. [S] 
Telluric blending, important in this part of the spectrum, is eas- 
ily identified by comparing the non-velocity-corrected spectra 
of both stars, as shown in Fig. [6] A reliable nitrogen abun- 
dance is even more difficult to obtain, since the only N line 
that is clearly present in our spectrum is weakly blended by Si I 
and S I lines. As a consequence, we were only able to deduce 
upper limits for nitrogen, being log e(N) < 6.6 for DS Aqr and 
log e(N) < 6.5 for V453 Oph. The derivation of these upper lim- 
its is shown in the synthesis in Fig|5] The oxygen abundance is 
primarily based on the strong (non-LTE sensitive) 7774 A mul- 
tiplet. For both stars this resulted in a slight underabundance 
of C and an overabundance of N and O, as reported in Table 
|5] This N overabundance is probably a consequence of the first 
dredge-up for both stars, while the O overa bundance is likely 
of initial composition (e.g. B oesgaard et alJ l999). 

Abundances are hard to obtain from sodium (Na) to sul- 
phur (S) due to the lack of lines. However most of the a- 
elements (Mg, Si, Ca and Ti) were detected, showing a mean 
of \a/ Fe\— 0.3 for both stars, which is expected for their metal- 
licity ( Lambertl l 19871 Idratton & Snedenlll99ll) . The Al abun- 
dance found for V453 0ph ([Al/Fe]= -0.59) differs signifi- 
cantly from this, but it is in agreement wi th the trend of the 
Al ab undance seen in metal-poor stars (e.g. lGratton & Snedenl 
1988). The intermediate mass elements from V to Zn are eas- 
ier to determine. These elements follow the Fe abundance 
more closely except for Cr and Mn. While for Mn the under- 
abund ance in both objects is expected (e.g. lGratton & Snedenl 
Il99lh . the rather large underabundance of Cr in DS Aqr is not. 




log e(N) = -« (and [S/Fe] = +0.2) 

DS Aqr: log e(N)=6.6, V453 Oph: log e(N)-6.5 

(and [S/Fe] = 0) 

j I , , , , I , , , I , i 

8675 8680 8685 

Wavelength (A) 

Fig. 5. Derivation of the upper limit for the nitrogen abundance, 
using spectrum synthesis of the N I triplet at 8680 A. Two syn- 
theses for each star were made: one without any nitrogen, to 
detect possible blends (dotted line), and the other one with the 
proposed upper limit abundance for nitrogen (thick continuous 
line). From the first synthesis, it turned out that the strongest 
N I line at 8680.28 A is not reliable to deduce an N abundance, 
due to weak Si I and S I blends. Therefore, the upper limits 
obtained in the second synthesis focus on the weaker line at 
8683.40 A, being log e(N) < 6.6 for DS Aqr and log e(N) < 6.5 
for V453 Oph . Logfgf) va lues for the three N I lines were taken 
from lHibbert et all Jl99ll) . 

A specific procedure was used to find all available lines for 
the s-process elements (Y, Zr, Ba, La, Ce, Nd and Sm). All lines 
in the VALD database were used to estimate their equivalent 
width for the assumed stellar parameters. All lines were then 
measured in the spectra according to their estimated strength 
up to the detection limit (EW~ 5 mA). Using this method, we 
maximised the number of lines per species and avoided mea- 
suring blends. Some of the s-process lines used in the analy- 
sis are shown in Fig. Only for Ba was the abundance es- 
timated using fewer than five lines. The abundance of Ce II, 
obtained using the a tomic data from the D.R.E.A.M. database 
dPalmeriet alfc OOO) is also reported. Using this data instead of 
the VALD data, the abundance of Ce remains fixed for DS Aqr, 
while it increases for V453 Oph. 

Table|5]gives the abundances of the light (Is) and heavy (hs) 
s-process elements. The light s-elements have a neutron num- 
ber close to the magic number 50 (Y and Zr), while the heavy 
s-elements cluster around the magic neutron number 82 (Ba, 
La, Ce, Pr, Nd and Sm). The indices indicate a mild overabun- 
dance of both Is as hs elements for V453 Oph, while for DS Aqr 
the s-process elements scale with Fe. 

The Eu abundance was determined for both objects. 
However, Eu lines are known to have a strong hyperfine split- 
ting (hfs). A reliable Eu abundance can therefore only be de- 
rived when this effect is taken into account. The Eu abundance 
of DSAqr is based on two EuII lines: the resonance line at 
4129.72A (W,i=94mA) and the line at 6645.06A (14mA); 
the Eu abundance of V453 Oph is derived solely from the res- 
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Synthesis of four CI lines 
DS Aqr: log e(C) = 6.80 
V453 Oph: log e(C) = 6.08 



9060 



9070 9080 
Wavelength (A) 



9090 



Fig. 6. Synthesis of the four C I lines used in the abundance analysis. These four lines turned out to be the only ones from which 
a reliable carbon abundance could be derived. The observed spectra are shown with a thin line; the syntheses with a thick line. 
The spectra are not velocity corrected, providing that the telluric lines are shown at their rest wavelength, and that they can be 
easily recognised. The macroturbulent veloc ity used in the synth eses is 8 km s for DS Aqr and 10 km s _I for V453 Oph. Log(g/) 
values for these four lines were taken from Hibbe rtet alJ Jl993h . 



onance line (24 mA). A synthesis of these lines including their 
hyperfine components showed that the weak lines are only af- 
fected by hfs in their profile, and not in their equivalent width. 
For the resonance line in DS Aqr, however, we found a small 
but significant effect: the difference in abundance between a hfs 
and a non-hfs treatment of this line is -0.06 dex. H yperfine A 
and B constants for this synthesis were taken from lLawler e t al. 
(2001b) and the solar isotopic composition of Eu was applied. 
The Eu abundance reported in Table [5] is the hfs-corrected 
abundance. A comparison between the Eu and Dy abundance 
give s [Dv/EvjI = 0-04 fo r both stars. This value is in agreement 
with Burris et alj j2000) where a [Dy/Eu] abundance very close 
to zero was found for all metallicities (-3.0 < [Fe/H] < -1.0). 

4.2. Error analysis 

Considering the absolute abundance determination (log(e)), the 
line-to-line scatter (<xi t i) represents the variance of the obtained 
values. Hence, if more than about five lines were used in the 



abundance determination, this scatter is representative of the 
non-systematic errors on the abundances. As indicated in Table 
13 the typical line-to-line scatter is 0.1 dex. 

Systematic and consequently more important errors on the 
derived abundances are mainly due: uncertainties on the deter- 
mination of atmospheric parameters, systematic errors in the 
log(gf) values and the presence of non-LTE effects. To deter- 
mine how sensitive the abundances are to the stellar parameters, 
the abundances were recalculated changing T e fr, log(g) and ( t 
with the assessed error (see Sect. 13. 3> . From this test, it appears 
that the abundance determination is most sensitive to the effec- 
tive temperature. The total error due to uncertainties in stellar 
parameters is given by adding all independent contributions: 

The systematic errors due to the possibly biased \og(gf) 
values and due to non-LTE effects are more difficult to deter- 
mine. Because of the similar effective temperature and gravity 
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Fig. 7. The velocity corrected spectra of the programme stars showing some of the used s-process lines. For V453 0ph, the 
s-process lines are almost the only visible lines in the spectrum. 



of both stars, a large reduction of these errors is guaranteed in 
a relative abundance analysis (see Sect. 15. 31 . 

In Table |5j the total error on the mean (cr tot ) of the [el/Fe] 
values is shown. This error is determined by adding the errors 
on the mean due: uncertainties in the model atmosphere (cr mo( j), 
the line-to-line scatter on the absolute abundances (defaulted to 
criti = 0.15 if only one line was used) and the error on the Fe 
abundance (cr Fe ). The employed formula is: 



0"tot 



(^) 2 + (cr mod )2 + (^p 



:) 2 



(2) 



5. Discussion of chemical analysis 

5.1. Comparison with the literature 

The chemical patterns of both stars were previously analysed 
in the literature. For DSAqr, two analyses were performed 
JGiridhar etalJll998l E"000). with the latter one being more 
elaborat e, while for V453Q ph only one analysis was per- 
formed JGiridhar et al.ll998l) . 

The analysis in the literature for V453 Oph and the one pre- 
sented here resulted in the same metallicity and very similar 



[el/Fe] values. Only for [C/Fe] and [O/Fe], being respectively 
0.5 dex lower and 0.8 dex higher in this analysis, there are dis- 
crepancies. Neither previous abundances are compatible with 
our spectra. For the other elements reported in both analyses, 
the mean difference relative to Fe is only 0. 12 dex. The high ra- 
dial velocity, indicating membership of pop. II stars, strength- 
ens the very low metallicity of the star V453 Oph. 



The metallicity reported in Giridh ar" et al.l (|2000) for 
DSAqr, [Fe/H]= -1.1, differs significantly from the one re- 
ported here. However our analysis uses a much larger number 
of lines of Fe I and Fe II , so we are confident in the established 
value. For the fractional abundances, however, the [el/Fe] val- 
ues in Giridhar et al. (2000) differ on average by only 0.06 dex 
from those reported here. Since we focused on s-process lines, 
we use many more available lines to derive quantitative data. 

We conclude that our spectra are clearly not compatible 
with some absolute abundance values found in the literature. 
Our S/N is higher and our abundance analysis covers many 
more elements, in particular more lines for s-process elements. 
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Table 5. Chemical composition of DS Aqr and V453 Oph as derived from the homogeneous abundance analysis. For every ion 
the abundance is listed absolute and relative to iron. D.R.E. A.M. means that the abundance of Ce II is determined using atomic 
data obtained from the D.R.E.A.M. database. The reported Eu abundanc e is corrected for hyperfine splitting. Also listed is the 
solar abunda nce (log(fr^)), taken from the same references as in iRevniers et alJ J2004) except for Dy, where we used the solar 
abundance of iBiemont & Lowd ([1993). Although more recent data exist, the values that were used, ensure that the employed 
logCgf) values are consistent with the solar abundances. The line-to-line scatter (cn t i) is reported if more than one line was used, 
together with the number of lines used to determine the abundance (N) and their mean equivalent width (W). Also reported is 
the error on the abundances relative to iron (cr tot ), determined as discussed in Sect. 14.21 This Table ends with the light and heavy 
s-process indices (respectively [ls/Fe] and [hs/Fe]) for the homogeneous analysis. The [ls/Fe] index is the mean of the Y and Zr 
abundance and the [hs/Fe] is computed using the mean of the Ba, La and Nd abundance. Clearly V453 Oph shows a significant 
overabundance for the s-process elements, while DS Aqr does not. 
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5.2. Depletion process 



As expected for RVC variables, the depletion process was not 
efficient in the investigated objects. This is clearly demon- 
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strated in Fig. [8] where the abund ances are plotted versus 
the condensation temperature jLodders & Feglevl 19881) . A star 
subjected to the depletion process shows a lowered abundance 
of elements with a higher condensation temperature. This is 
clearly not observed in either variable. Therefore, the derived 
abundances represent the initial composition of the stars, pos- 
sibly altered by their internal evolution and/or accretion of en- 
riched material from a companion. 



Table 6. Results from the relative analysis of DS Aqr and 
V453 Oph only using the lines that were observed in both spec- 
tra. In the column labelled A[el/Fe], the difference in abun- 
dance relative to Fe for V453 Oph minus the one for DS Aqr 
is given. The other notations are similar as in Table H] This 
analysis shows clearly an overabundance of s-process elements 
in V453 0ph because possible bias effects due to uncertain 
log(gf) values are eliminated in a relative study. 

Relative 
Analysis 
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Fig. 8. The derived abundances versus the condensation tem- 
perature for DS Aqr (top) and V453 Oph (bottom). These fig- 
ures indicate that no depletion process occurred. 



5.3. s-process overabundances 

The s-process indices reported in Table|5]in combination with 
the abundances of the individual s-process elements indicate 
that DS Aqr shows the expected abun dances for a photosp here 
unaltered by internal evolution (e.g. lFrancois et alJl2003l and 
references therein), while for V453 Oph the overabundance of 
heavy elements indicates mixing of nucleosynthetic elements 
in the photosphere. 

To constrain possible systematic uncertainties in the re- 
ported abundances, a strictly relative analysis was performed, 
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which uses only spectral lines detected in both objects. The 
results of this analysis are listed in Table [5] By comparing 
V453 Oph with the results of DS Aqr, we assume that the lat- 
ter is a good tracer of the mean Galactic trend. The individual 
relative abundances given in Table[6](A[el/Fe]) are thus to be 
interpreted as the abundance characteristics of V453 Oph with 
respect to unevolved objects at similar metallicity that follow 
the mean Galactic trend. This analysis provides a solid basis 
for the proposition that V453 Oph is a metal deficient pop II 
RV Tauri star with a significant enrichment in heavy elements. 

The errors on the mean of the relative abundances (cr tot ) 
were estimated assuming that uncertainties due to the model 
atmospheres cancel out because both stars have similar atmo- 
spheric parameters. Thus we used a formula similar to Eq. |2 
excluding the uncertainty on the model atmospheres, but in- 
cluding the uncertainty on the Fe abundance of both stars. The 
resulting errors on the relative abundances are shown in Table 

E 
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V453 Oph 



[Fe/H] 

Fig. 9. In a narrow metallicity interval, the literature [hs/ls] val- 
ues and those for V453 Oph are plotted. We refer for the def- 
inition of the Is and hs indexes to the text. The CH stars with 
the expected Pb overabundance are depicted with filled sym- 
bols, while the others ar e represented using open symbols. The 
circles ar e taken from I Van Eck et all J2003I) . the upper trian- 
gles fro m lAoki et al.1 fcOOlh and the squares from Ao ki et alJ 
il2002l) . The star V453 Oph is shown as an asterisk and labelled 
in the graph. For this star, no Pb abundance was determined. 



5.4. Comparison with CH-giants 

We performed a comparison of the V453 Oph s-enrichment 
with those reported for CH stars of similar metallicity. CH-stars 
are cool C-rich objects of low metallicity showing enhanced 
s-process lines. A distinction is made between early-type and 
late-type CH-stars: Early-type CH-sta rs are probably binaries 
which are extrinsically enriched (e.g. McClure & Woodsworth 
1990), while late-type CH-stars combine bright absolute mag- 
nitudes and high 12 C/ 1 3 C ratios and are ther efore believed to be 
intrinsically enriched (Van Eck et al. 2003). 

Standard s-process scenarios predict the production of 
mainly Pb for the low-metallicity CH-stars (Pb/Ba ratios of 7 
are predicted for [Fe/H]= -1.3 bv lGorielv & MowlavilhoOOl) . 
Some of the studied CH-st ars show this extreme Pb enrichment, 
while others do not (e.g. I Van Eck et alJl2003l and references 
therein). A comparison of the reported s-process enrichments 
is shown in Fig. [9] In this figure, the [hs/ls] value is plotted 
as a function of metallicity. The Is value is taken as the abun- 
dance of Zr and the hs value as the mean of the abundances 
of La, Ce and Nd. Only stars with all these abundances deter- 
mined were taken into account. In Fig.[9]a large spread in [hs/ls] 
is observed for very similar metallicities. Although V453 Oph 
is situated in a less dense region of this plot, some CH-giants 
have similar [hs/ls] indices. Extremely Pb-enriched as well as 
not Pb-enriched objects are found with a [hs/ls] value similar 
to V453 Oph. Note that the determination of the Pb abundance 
in the hot V453 Oph star is difficult and only an upper limit of 
[Pb/Fe] < 2 could be determined. 



5.5. Other RVTauri stars with s-process enrichment 

To our knowledge, up to now only one other RV Tauri object 
with a clear s-process en richment has been detected: VI in w 
Cen dGonzalez & Wallersteir]ll994 . This RV Tau star with a 
spect roscopic class on the B-C boundary dWallerstein & Coxl 
1984) also shows mild s-process overabundances ([ls/Fe] = 0.6 
and [hs/Fe] = 0.5). It is unclear whether this enrichment is in- 
trinsic, or if it is due to the well-studied bu t poorly un derstood 
self-enrichment of o> Cen ("e.g. lTsuiimoto & Shigevamall2003l 
and references therein). The fact that t he less evolved variable s 
V48 and V29, also studied in lGonzalez & Wallersteir] |l994), 
display the same s-process overabundances as VI, is in favour 
of the self-enrichment scenario. VI, however, also shows a 
small, but significant, overabundance in carbon compared to 
these two other variables (A[C/Fe]=0.9 and 0.5 for V48 and 
V29 respectively). Thus it appears that it is not only possible to 
enhance solely the s-process elements (V453 Oph) in RV Tauri 
variables, it is also possible to enhance only the C abundance 
(u) Cen VI) without strongly enriching the star with s-process 
elements. 



6. V453 Oph, initial, extrinsic or intrinsic heavy 
element overabundances ? 

The most important abundance characteristic of V453 Oph is 
that it is a genuine low-metallicity star with [Fe/H] =-2.2 show- 
ing a significant heavy element overabundance, but without any 
C-enrichment. Moreover, the N-abundance is too low to ex- 
plain the C deficiency by CN cycling. The luminosity and at- 
mospheric parameters of V453 Oph are compatible with an ob- 
ject of low initial mass (typically M =t 0.8 M ) in its post-AGB 
phase of evolution. 

We explored several possible scenarios to explain the ob- 
served abundances: special chemical characteristics of the 
parental cloud; enrichment by a binary companion and intrinsic 
s-process enrichment by dredge-up. 

6.1. Parental cloud anomalies 

At low metallicity, the Galaxy was not well mixed so there is 
a large star-to-star scatter in the initial abundance ratios of un- 
evolved stars of low iron abundance ([Fe/H] < -2). The chemi- 
cal composition of low-metallicity stars may reflect strong ini- 
tial heavy element enrichment of the parental cloud. To test 
whether this may have been the case for V453 Oph, we scanned 
the recent literature for good quality data of individual non- 
evolved stars. We compared both the overabundances of indi- 
vidual neutron capture elements as well as the observed distri- 
bution, since the latter is a good tracer of the nucleosynthesis 
production site. 

Massive stars are believed to be major contributors of neu- 
tron capture elements in the Galaxy, with a most probable en- 
richment in r-process nuclei produced during the supernova ex- 
plosion, but also of the so-called weak s-process component 
responsible for the production of s-nuclei up to A 90 during 
core He-burning. Low-metallicity stars polluted by supernova- 
enriched material are therefore likely to show, in addition to 
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La (Simmerer) 
Eu (Simmerer) 
V453 Oph 



[Fe 



2.0 
ll/H] 



Table 7. Comparison between the enrichments found in 
V453 0ph and the abundan ces found in unev olved low- 
metallicity stars compiled by iBurris et"a"fl J2000l) . The mean 
and standard deviation of the values obtained for the 18 stars 
in the metallicity range -2.4 < [Fe/H] < -2.0 is given. Not all 
data are available for all stars, so the number of stars used in 
the compilation is given in the last column. 
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Fig. 10. A comparison of the La and Eu abundance of 
V453Qph with the va lues found for non-enriched objects in 
ISimmerer et al.l (Eo04). Although the excess is small, the La 
abundance of V453 Oph is above the values observed for non- 
enriched objects. 

an s-process enrichment in elements up to Zr, a strong r- 
process signature. In recent years, it became clear that this pic- 
ture is oversimplified and recent observations of a large sam- 
ple of low-metallicity stars show evidence of the main and 
strong components of the s-process in individ ual unevo lved 
stars with metallicities as low as [Fe/H] =-3 dSimmerer et alJ 
120041) . These authors focus on the study of La and Eu abun- 
dances of about 160 giant and dwarf stars in the metallicity 
range -3<[Fe/H]<+0.3. In Fig[TO|we compare the abund ances 
of V453 Oph with the sample of Sim merer et alJ (|2004). The 
La enrichment is above the level observed in objects of similar 
metallicity, but the excess is mild. 

Detailed chemical studies of very low-metallicity un- 
evolved stars tend to focus on the few objects which are known 
to show this clear dominant r-process signature. The most fa- 
mous example in this respect is t he ultra-metal-deficient star 
CS 22892-052 JSneden et al.l2003[ and references therein). We 
compared the heavy element distribution of V453 0ph with 
the well studied r-p rocess enriched halo giant BD+17".3248 
dCowan et al.l E002). Both stars have the same metallicity and 
show similar overabundances of some individual heavy ele- 
ments. The difference between both distributions is most ap- 
parent in Eu and Dy. When both distributions are normalised 
to Zr, both the Eu and Dy abundances of BD+17".3248 are 0.6 
dex larger than what is observed in V453 Oph. 

In Table we compare the V453 0ph overabundances to 
those found in non-evolve d objects of simila r metallicities 
from the extensive study of IBurris et al.1 J2000l) . At metallic- 
ity [Fe/H]=-2.2, the Y abundance of V453 Oph ( [Y/Fe]=+0.3 
and [Y/Fe]= +0.5 in DSAqr) is clearly above the mean value 
of -0.2 from lBurris et alJ ( | 2000) sample. Outliers to this mean 
Galactic trend exist, however, but the highest observed value 
in this sample of 10 stars is about [Y/Fe]=+0.0 in the metal- 
licity range of V453 Oph. The scatter in the [Zr/Fe] ratio at 



metallicity around -2.2 is similar, but here again V453 Oph is 
significantly more Zr-rich than the other objects. The La over- 
abundance with respect to the Galactic mean is confirmed in 
this comparison. For the heavy elements Nd, Eu and Dy, the 
general picture is not as clear. In V453 Oph, they tend to be as 
abundant as in the IBurris et"afl d2000l) sample, although some 
stars do exhibit similar enrichments. 

In conclusion, although the [Eu/Fe] and [Dy/Eu] abun- 
dances in both V453 Oph and DS Aqr indicate an initial enrich- 
ment of r-process nature, the relatively large Y, Zr, Ba and La 
overabundances of V453 Op h with respect to DS Aqr and the 
low-metallicity stars in the Burri s et alJ (J2QO0) sample clearly 
points to an additional s-process enrichment in V453 Oph. If 
the actual chemical distribution of V453 Oph reflects the initial 
composition of the parental cloud, this cloud must have been 
subject to some s-process enrichment capable of producing el- 
ements up to the Ba peak. However, up to now, only AGB stars 
are believed to be responsible for this corresponding main com- 
ponent of the s-process, massive stars being capable of produc- 
ing elements up to Zr only. 

6.2. Extrinsic enrichment 

The enrichment we see in V453 Oph may come from the pol- 
lution of a putative companion similar to the extrinsically en- 
riched CH stars discussed in Section 5.4. Binarity among dusty 
RV Tauri stars may be widespread dVan Winckel et al1 [l999b; 
De Ruvter et alJ2005l) and the arguments are either direct (with 
direct radial velocity measurements), or indirect (SED charac- 
teristics pointing to the presence of a dusty disc instead of an 
outflow). Note that only the dusty RV Tauri stars are thought 
to be binaries in which the companion is probably not a white 
dwarf but an unevolved main sequence star. No binary RV Tauri 
star is definitely identified with a White Dwarf companion. 

Unfortunately, we do not have detailed radial velocity mon- 
itoring data for V453 Oph. In V453 Oph, no dust-excess is ob- 
served so the above mentioned indirect evidence for binarity is 
lacking. Also, the intrinsic luminosity determined by the P-L 
relation is large, pointing to a genuine post-AGB nature. Note 
that extrinsically enriched objects are rare : Giant barium stars 
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(extrinsic ally enriched Giants with white dw arf companion s) 
represent only 1 % of the giant G-K stars (e.g. Jorissen 2004). 

We conclude that, although we cannot exclude the possibil- 
ity that V453 Oph belongs to a binary system which has been 
polluted by an AGB star that is now a White Dwarf companion, 
this hypothesis is not very likely. Moreover, the C deficiency 
associated with the s-process enrichment remains difficult to 
explain by standard nucleosynthesis models (see Sect. 7). This 
scenario presents however the advantage that the White Dwarf 
progenitor responsible for the s-process nucleosynthesis can a 
priori be of any initial mass (below some 8 M ). 

6.3. Intrinsic enrichment 

The most natural scenario for a post- AGB star to become over- 
abundant in s-process elements is that the object was subject 
to internal nucleosynthesis followed by dredge-up episodes 
during its previous AGB phase. The most favoured s-process 
model is associated with the partial mixing of protons (PMP) 
into the radiative C-rich layers at the time of the 3D UP 
jGallino et alll998UGorielv & Mowlavl2000l) . There is, how- 
ever, no general agreement on what physical mechanism can be 
responsible for the injection of protons into the 12 C-rich layers. 
Furthermore, in this scenario, the s-process surface enrichment 
is bound to be associated with a stronger C pollution. 

With a metallicity of -2.2 and a luminosity L = 2400 ± 
1600 L Q , this post-AGB star must have evolved from a low- 
mass main sequence star (typically M ^ 0.8 M ). For such 
stars, it remains unknown if any PMP may take place in the 
C-rich region, and even if any third dredge-up can intrinsically 
enrich the stellar surface. The results of our abundance analysis 
are consequently difficult to reconcile with this evolutionary 
path as well. 



7. S-process models and observations 

As explained in the previous section, none of the standard sce- 
narios can a priori explain the present observations. It is clear 
that the abundances of V453 Oph challenge our current un- 
derstanding of the s-process. On the one hand, s-process nu- 
cleosynthesis during core He-burning in massive stars can not 
account for a significant production of s-elements up to the Ba 
peak. On the other hand, any AGB s-process models within 
the PMP scenario would lead to a concomitant and significant 
production of C whic h is not observed. More specifically, as 
described in detail in JDeroo et alJl2005h , the PMP model in 
the intrinsic enrichment scenario could in principle lead to an 
abundance distribution close to the observed one, but would 
simultaneously pollute the stellar envelope with a C overabun- 
dance of the order of [C/Fe]=1.5. The major specificity of this 
PMP s-process is the prediction of a large Pb overabundance 
that cannot be confirmed by the present observations (the up- 
per limit established from the observations is [Pb/Fe]< 2.0). 

As shown by ( Der oo et alJl20 05). there is however one s- 
process model that could explain to some extent at least part 
of the observations. This is the so-called convective s-process 
model that can possibly develop at the base of hot thermal 





1 




0.8 




0.6 








0.4 






0.2 









-0.2 



o V453 0ph 
— Convective s-process in a 3Mo 

Bal 




40 



50 



60 

z 



70 



80 



Fig. 11. Comparison of the observed surface abundances (en- 
richment of V453 Oph with respect to DS Aqr) with the s- 
process abundances obtained within the convective s-process 
model in a 3M Z=0.0001 model star. 



pulses in massive AGB stars. Assuming again the extrinsic en- 
richment scenario to prevail, we can imagine that a massive 
AGB companion star could have polluted V453 Oph by stellar 
winds. In massive AGB stars, the temperature at the base of the 
thermal pulse is hot enough for the 22 Ne(a,n) 25 Mg to become 
active and be responsible for some s-process nucleosynthesis. 
The major advantage of this scenario in comparison with the 
PMP model is that the s-process enrichment in the pulse can 
reach values similar to or even greater than the C one. This is 
due to the large mass extent on which s-process seed nuclei are 
irradiated by neutrons (M pu i se ~ 1O~ 3 M ) compared with the 
lO~ 5 M resulting from the PMP models. The resulting surface 
overabundances are shown in Fig. ^2 Such a dilution gives rise 
to the corresponding overabundances [C/Fe]=0.1, [N/Fe]=0.6 
and [Mg/Fe]=0.4, and no Pb overproduction contrary to the 
PMP nucleosynthesis. There is however a drawback to this 
scenario: to reduce the C pollution, the star should not have 
endured the third dredge-up episodes. It remains unclear what 
star may or may not be subject to the third dredge-up, different 
convective models with or without extra mixing mechanisms 
leading to significantly different predictions. The presence of a 
binary companion is an additional source of complexity. 

Only in this very specific model and if no 3DUP occurs 
to enrich the envelope in carbon, the surface enrichment of s- 
elements can be reconciled with a low C-content. It is clear 
that, although the comparison of the calculated and observed 
abundances is quite satisfactory, the 3M Z=0.0001 model star 
can only have been the putative companion to V453 Oph. Such 
a massive object of low metallicity must have evolved into a 
white dwarf by now. 

8. Conclusion 

In this paper, an abundance analysis is presented of two very 
similar RV Tauri stars of class C, V453 Oph and DS Aqr. Both 
objects are genuine post-AGB stars, although both have no de- 
tected IR excess due to circumstellar dust. These stars differ 
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significantly when we look at their photospheric abundance 
pattern. While V453 Oph shows an excess of light and heavy 
s-process elements, DS Aqr shows neither. This conclusion is 
clearly demonstrated by the relative analysis that was per- 
formed between both objects. 

The abundances derived for V453 Oph are of particular nu- 
cleosynthesis interest, since for the first time, the observed 
overabundance of a s-process signature is not accompanied 
by a simultaneous C-enrichment. Moreover, the N abundance 
upper limit excludes possible CN cycle destruction of the 
dredged-up carbon. The most favoured AGB s-process model 
based on the partial mixing of protons at the time of the third 
dredge-up is inevitably accompanied by a larger surface enrich- 
ment in C and therefore fails to account for the present obser- 
vation. 

We explored three evolution scenarios to explain the abun- 
dance distribution of V453 Oph. All have specific shortcom- 
ings. The first possibility is that the parental cloud was initially 
polluted by heavy elements. The abundances observed in non- 
evolved objects of the same metallicity ([Fe/H]~ -2.2) show 
quite some spread, illustrating that the Galaxy was not well 
mixed at these epochs. Although some of these objects do ex- 
hibit significant heavy element overabundances, a comparison 
with the abundances observed in V453 Oph indicate that its 
parental cloud must have been quite extremely polluted in r- 
but also s-elements. 

Alternatively, the star may have been enriched during its 
evolution either by mass transfer or internally. If the star is in- 
deed intrinsically enriched, this would suggest that low-mass 
objects of low metallicity can be subject to s-process nu- 
cleosynthesis without simultaneously polluting their envelope 
with C. At the present time, it is impossible to reconcile such 
conditions with our present understanding of AGB stars. If the 
star belongs to a binary system, the companion may have pol- 
luted the actual post-AGB star but here again the same C de- 
ficiency is difficult to reconcile with the traditional s-process 
scenarios. Moreover, there is still no evidence of a WD com- 
panion. Nevertheless, we explored the possibility for a 3M 
Z=0.0001 model star to produce s-elements within the ther- 
mal pulses only and conclude that such a convective s-process 
model could to some extent reconcile the observed s-process 
enrichment without significant C-enrichment. Clearly, more in- 
depth theoretical and observational studies are needed to fur- 
ther understand if such a scenario could indeed be responsible 
for the extraordinary chemical composition of the V453 Oph 
surface. 
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